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Based on the finite time thermodynamics theory, the entransy theory and the entropy theory, the Stirling cycles under different 
conditions are analyzed and optimized with the maximum output power as the target in this paper. The applicability of entran¬ 
sy loss (EL), entransy dissipation (ED), entropy generation (EG), entropy generation number (EGN) and modified entropy 
generation number (MEGN) to the system optimization is investigated. The results show that the maximum EL rate corre¬ 
sponds to the maximum power output of the cycle working under the infinite heat reservoirs whose temperatures are prescribed, 
while the minimum EG rate and the extremum ED rate do not. For the Stirling cycle working under the finite heat reservoirs 
provided by the hot and cold streams whose inlet temperatures and the heat capacity flow rates are prescribed, the maximum 
EL rate, the minimum EG rate, the minimum EGN and the minimum MEGN all correspond to the maximum power output, but 
the extremum ED rate does not. When the heat capacity flow rate of the hot stream increases, the power output, the EL rate, the 
EG rate and the ED rate increase monotonously, while the EGN and the MEGN decrease first and then increase. The EL has 
best consistency in the power output optimizations of the Stirling cycles discussed in this paper. 
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1 Introduction 

Under the background of the global energy shortage and the 
environment pollution, the exploration for the efficient and 
clean energy utilization becomes one of the focuses of the 
attention now. The Stirling engine, which is also called the 
external combustion hot-gas engine, has great potential for 
the future applications due to its high efficiency and good 
adaptive capacity for various kinds of energy and fuels [1]. 
Therefore, the optimization of the Stirling cycle is of great 
significance. 

Curzon and Ahlborn [2] first proposed the concept of fi¬ 
nite time thermodynamics (FTT) in 1975. As an extension 
of the traditional irreversible thermodynamics, the essence 
of the finite thermodynamics studies is to reduce the irre- 


^Corresponding author (email: liangxg@tsinghua.edu.cn) 


versibility and to optimize the thermal performances of 
practical systems by introducing heat transfer theory into 
thermodynamic analyses [3]. The FTT theory has been 
widely used by many researchers to analyze and optimize 
the energy systems, including heat engines, refrigerators, 
heat pumps, etc., under the constraints of finite heat transfer 
rate or finite size [4-25], promoting its rapid development 
[3, 26, 27]. As a typical model of the FTT analyses, the 
Stirling engine has received great attention. Li et al. [28] 
established the Stirling engine model with finite heat trans¬ 
fer rate, heat leak and finite regeneration time, and investi¬ 
gated the power and efficiency optimizations of the solar- 
powered Stirling heat engine by the FTT theory. They found 
that the absorber temperature and the regenerator effective¬ 
ness had significant effects on the performance of the Stir¬ 
ling engines. More FTT researches can be found in refs. [1, 
29-34], etc. 


© Science China Press and Springer-Verlag Berlin Heidelberg 2012 


tech.scichina.com www.springerlink.com 




229 


Zhou B, et al. Sci China Tech Sci January (2013) Vol.56 N0.1 


Entropy generation (EG) is usually regarded as the 
measure of the ability loss of doing work in thermodynamic 
cycles. Smaller EG indicates less loss of available energy. 
Therefore, the minimum EG is related to the optimal per¬ 
formance in thermodynamic cycle optimizations [35, 36]. 
Bejan [36] developed the EG minimization theory that takes 
the minimum EG as the optimization objective in consider¬ 
ing the effect of finite temperature difference heat transfer 
on the total EG of the thermal system. Based on this theory, 
the minimum EG principle has been widely applied to the 
optimizations of thermal processes, such as heat-work con¬ 
version processes [11, 37-39]. However, arguments arose in 
its application. Some researches show that the minimum EG 
does not always correspond to the optimal performance of 
the thermal systems [40, 41]. The applicability of the mini¬ 
mum EG principle needs further discussion. 

The entransy theory is newly developed in recent years 
and is applied to heat transfer optimizations. The concept of 
entransy is proposed by Guo et al. [42] from the analogy 
between the heat and electric conduction. Based on this 
concept, Guo et al. [42] developed the extremum ED prin¬ 
ciple and the minimum entransy-dissipation-based thermal 
resistance principle. These principles are used to optimize 
heat conduction [43, 44], heat convection [45-47], thermal 
radiation [48, 49] and the design of the heat exchangers 
[50, 51]. For the concept of entransy dissipation (ED), re¬ 
searches show that it is not appropriate for the optimization 
of heat-work conversion processes [52, 53]. However, ef¬ 
forts are made to extend the application of the entransy the¬ 
ory to the optimization of the heat-work conversion. Cheng 
et al. [41, 54] proposed the concept of entransy loss (EL), 
which is the sum of ED due to the irreversible heat transfer 
and the entransy variation due to the work output. EL is the 
entransy that is used in the heat-work conversion processes. 
The analyses of the air standard cycle [41] and the air con¬ 
ditioning system [54] show that the increase in the EL rate 
always leads to the increase in the power output. Compared 
with the optimization based on the concept of entropy, there 
is a better consistency between the variations of EL and the 
power output. However, there are not many researches re¬ 
lated to EL, and no reports about the optimization of the 
thermodynamic cycles by combining FTT and EL currently. 

Based on the FTT theory, the entransy and entropy theo¬ 
ries are applied to the analyses of the Stirling cycles under 
different conditions in the present paper. The maximum 
power output is the objective. The applicability of the con¬ 
cepts of EL, ED, EG, etc., in the FTT optimization is dis¬ 
cussed. 

2 Model of the Stirling cycle 

The typical simplified flow chart of the Stirling engine and 
its T-s diagram are shown in Figure 1. The Stirling engine is 
mainly composed of the compression space, the expansion 


space and the regenerator, etc. There are four processes in 
the Stirling cycle as shown in Figure 1: (1) a-b, the isother¬ 
mal process in which the working fluid absorbs heat and 
expands for doing work; (2) b-c, the isochoric process in 
which the working fluid keeps its volume and transfers heat 
to the regenerator; (3) c-d, the isothermal process in which 
the working fluid is compressed and releases heat to keep its 
temperature; (4) d-a, the isochoric process in which the 
working fluid keeps its volume and absorbs heat from the 
regenerator. For ideal regeneration, the heat absorbed dur¬ 
ing process d-a equals the heat released during process b-c. 
However, under the limits of finite regenerator size and the 
finite regenerate rate in practical applications, the ideal re¬ 
generation process cannot be realized. Therefore, the effect 
of the regenerate loss AQr is considered in the present pa¬ 
per. 

For the Stirling cycle shown in Figure 1, assume that the 
working fluid is the ideal gas whose mass is n mol, and the 
working fluid temperatures at a-b and c-d processes are 7) 
and r 2 , respectively. Then, the heat absorbed in process a-b 
and the heat released in process c-d by the working fluid are 
[35] 


Q l = nR'I\ In A, 

0) 

Q , = nRT 2 In A, 

(2) 


where R (J (mol K) _1 ) is the universal gas constant, and 
i=14/Va is the compression ratio. Assume that the effec¬ 
tiveness of the regenerator is r] m , the regenerated heat is [28] 



(a) 



(b) 


Figure 1 Flow chart and T-s diagram of the Stirling engine, (a) Typical 
simplified flow chart of the Stirling engine; (b) T-s diagram of the Stirling 
cycle. 
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Qr = nc v riJT l -T 2 ). (3) 

The regeneration loss due to the imperfection of the re¬ 
generation process is [28] 

AQ R =nc v (l-TjJ(T 1 -T 2 ). (4) 

The existence of the regeneration loss is equivalent to a 
heat transfer of AQr that transfers directly from the hot 
working fluid to the cold working fluid. Therefore, the total 
heat supplied by the hot reservoir <2h and that received by 
the cold reservoir Q L in a cycle are 

Qn = Gi + AG r = nRT x In A + nc v (1 - t] m )(7] - T 2 ), (5 ) 

Q l =Q 2 + AQ r = nRT 2 In/, + nc v (1 - rj m )('!] -T 2 ). (6) 

The following sections will show the analyses of the 
Stirling cycles with the maximum power output as the ob¬ 
jective based on the FTT theory. We will discuss two dif¬ 
ferent heat supplies, the infinite heat reservoir and finite 
heat reservoir, respectively. 

3 FTT optimizations under heat sources of in¬ 
finite heat reservoirs 

As an external combustion engine, the Stirling engine can 
have various kinds of heat sources [35]. When the heat res¬ 
ervoirs are provided by vapor, etc., the heat reservoirs can 
be considered to have infinite heat capacities whose tem¬ 
peratures are constant. The T-s diagram of this cycle is 
shown in Figure 2. 

Assume that the thermal conductances of the heat ex¬ 
changers on the hot and cold sides are (UA) H and (UA) L , 
respectively. The heat absorbed from the hot reservoir and 
that released to the cold reservoir by the working fluid dur¬ 
ing each cycle are 


Gh =(UA) h (T h -7j)r p 

(7) 

Ql - (UA) l (T 2 —T l )t 2 , 

(8) 


where ri and z 2 are the times of process a-b and process c-d. 



Figure 2 T-s diagram of the Stirling cycle under infinite heat reservoirs. 


respectively. The combination of eqs. (5)—(8) gives 

r _ nRT t In A + nc v (1 - r/ m )(7j - T 2 ) 
(UA) h (T h -7j) 

T _ nRT 2 In A + nc v (1 - ij m )(7j - T 2 ) 

(UA) l (T 2 -T l ) 


( 9 ) 

( 10 ) 


Considering that the regeneration processes are executed 
with finite heat transfer rate and assuming that the tempera¬ 
ture of the working fluid varies at a constant rate, the times 
of the regeneration processes b-c and a-d can be expressed 
as [28] 


1 

ii 

(ID 

z-4=(7]-7;)/M 2 , 

(12) 

where and M 2 are the constants related to the properties 

of the regenerator such as materials, structures, etc. 

Based on the analyses above, the cycle time of the Stir¬ 
ling cycle can be obtained 

r = r, + r 2 +r 3 +r 4 , 

(13) 

and the cycle output power is 


w = (Gh - Gl )/t = nR \nX{T x -T 2 )/t. 

(14) 


According to its definition, EL is the difference between 
the entransy flows that get into and out of the system [54], 
Dividing EL of each cycle by the cycle time, the EL rate of 
the cycle is 

G l0SS -(G H 7H-G L r L )/r. (15) 


EG and ED in a cycle are the respective sums of the EG 
and ED in all the processes. Therefore, the EG rate and the 
ED rate of the engine are 


.S' 


gen 


_'6 h Gh,G l G l , aq r 

l T t T h T l T 2 T 2 

= (Ql/t l -Q b /t h )/t, 



(16) 


d dls = [Qn (T h ~ 7]) + Q l (T 2 -TJ + AQ r (7j - T 2 )]/r 

= (Q h T h -Q l T l -Q 1 T 1 +Q 2 T 2 )/t. (17) 


It can be seen from eqs. (5)-(14) that the variations of 
the working temperatures 7j and T 2 will affect the output 
work and the cycle time at the same time when the temper¬ 
atures of the hot and cold streams 7 H and 7] are prescribed, 
and will consequently lead to the variation of the output 
power. For the objective of the maximum power output, the 
working temperatures 7) and T 2 need to be optimized. 

Assume that the working temperatures 7] and 73 and the 
heat reservoir temperatures T H and /] satisfy [34] 

TJT, ={TJT U )\ 


(18) 
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where r is the isochoric temperature ratio exponent. The 
objective of the present problem is to find the optimal 7j 
that leads to the maximum power output under arbitrarily 
prescribed r, and to investigate the applicability of the en- 
transy and entropy theories to the optimization. 

Let us consider a numerical example. Assume that r = 
0.6, T h = 1000 K, 7j,= 300 K, (UA) H = 300 W K 1 , ( UA) L = 
300 W K 1 , tj m = 0.8, c v = 20 J (mol Kf 1 , R = 8.314 
J (mol K) _1 , X = 3, n = 0.5 mol and M\= M 2 = 8xl0 4 K s -1 . 
The practical parameters are different for various engines due 
to the difference in the structures, working fluids, etc. As the 
present work discusses the applicability of the entransy and 
entropy to the optimization of the Stirling cycles, the above 
parameters are not originated from a real Stirling engine, but 
are selected in the general range of the practical engine pa¬ 
rameters, and are similar to values in refs. [28, 31]. 

The output power, the EL rate, the EG rate and the ED 
rate with different working temperature 7j are calculated 
based on the above parameters, and are shown in Figure 3. 
It can be seen that all the parameters increase first and then 
decrease with increasing 7j. The output power, the EL rate 
and the EG rate reach their maximum values when 7j = 817 
K. Therefore, the maximum EL rate corresponds to the 
maximum output power for this problem. The minimum EG 
rate does not correspond to the maximum output power, 
but the maximum EG rate does. The ED rate reaches its 
maximum value at T\ = 782 K, which means that the extre¬ 
mum ED rate does not correspond to the maximum output 
power. Thus, the ED rate is not suitable for the output pow¬ 
er optimization. 

The above results are analyzed below. The expression of 
the EL rate for this cycle can be changed into 

6 . _ Q H T H -Q L T L _ W(T H -T L )/ n+ WT L 

w loss 

T T 

W[T h -(1-v)T l ] ' 

n 

where // is the cycle efficiency. Similarly, the expressions of 
the EG rate and the ED rate also can be changed into 



7](K) 


Figure 3 Variations of the output power, the EL rate, the EG rate and the 
ED rate with T\. 


s 


gen 


w[q->7)r H -r L ; 


( 20 ) 


G dis = 


T H -(l->i)T L +riT 2 


W + 


ejr.-rj 


( 21 ) 


For the Stirling cycle working under infinite heat reser¬ 
voirs, the cycle efficiency can be calculated based on eqs. 
(5) and (6), 

>7 = (Gh-Gl)/2h 

nRT t In A- nRT\ In A 
nRT ] In A + nc v (\-7) m )(T t -T 2 ) 

= _ 1 ~TjT x _ 

l + [c y (l-TlJ(X-T 2 /T 1 )]/(R]nA.y (22) 


This equation and eq. (18) show that the cycle efficiency 
is constant when T H , T L and r are prescribed. Eqs. (19)—(21) 
show that the EL rate and the EG rate increase monoto¬ 
nously with increasing output power when the cycle effi¬ 
ciency is prescribed, while the ED rate has no direct relation 
to the output power under the same condition. 

From the viewpoint of entransy, there are three necessary 
processes in heat-work conversion, and the corresponding 
costs must be paid for each process. Firstly, heat is trans¬ 
ferred from the hot reservoir to the working fluid, and the 
corresponding cost is ED. Secondly, heat is converted to 
work by the thermodynamic cycle. Heat entransy is con¬ 
verted to work entransy in this process that leads to the de¬ 
crease of heat entransy. Finally, the used heat is transferred 
to the cold reservoir, and the cost of this process is also the 
corresponding ED. Overall, EL is the necessary costs of the 
systems for outputting work. Therefore, the EL rate can be 
applied to the optimization of the output power. However, 
ED does not describe all the processes of the heat-work 
conversion. That is why the ED rate cannot be used for the 
output power optimization. 

Bejan [36] introduced the concept of entropy generation 
number (EGN) when EG is used for the optimization of the 
thermal processes involving heat exchangers. EGN is 
expressed as 

N s = S gen /C mm , (23) 

where C m ;„ is the minimum heat capacity flow through the 
heat exchangers. The “EG paradox” is noted when EGN is 
used in the analyses of the counter flow heat exchangers 
[36]. To remove the paradox, Xu and Yang [55] defined the 
modified entropy generation number (MEGN), which is the 
production of the EG per unit heat transfer rate and the am¬ 
bient temperature. The expression is 

^ gs =U/2• ( 24 ) 

For the present problem, EGN cannot be defined because 
there is no heat capacity flows through the heat exchangers. 
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For MEGN, combining eqs. (20) and (24) yields, 

^ g , = r 0 [(i-» 7 )rH-r L ]/(r H r L ). (25) 

Eq. (25) shows that MEGN is a constant when T H , 7), T 0 
and r are prescribed. The output power has an extremum 
value with the variation of the working temperature 7) in 
such condition. So, MEGN cannot be used for the output 
power optimization in such prescribed conditions. 

The case in which r = 0.6 is discussed above. Further 
analysis of the above equations shows that the relationship 
remains the same for arbitrary values of r. The maximum 
EL rate and the maximum EG rate always correspond to the 
maximum output power with infinite heat reservoirs of con¬ 
stant temperature, while the ED rate and the MEGN have no 
direct relation with the output power. 

4 FTT optimizations under heat sources of fi¬ 
nite heat reservoirs 

The Stirling engine has great potential for the applications 
in the field of the waste heat recovery due to its high effi¬ 
ciency [35]. The used water, gas, etc., at certain temperature 
are the common waste heat sources in industry, and the heat 
can be transferred to the Stirling engines through heat ex¬ 
changers to produce power output. The flow chart and T-s 
diagram in such an application are shown in Figure 4. The 
hot stream at temperature T’e-m flows through the exchanger 
on hot side with heat capacity flow rate Ch, and the cold 
stream at temperature T L . in flows through the exchanger on 
cold side with heat capacity flow rate Cl- The used streams 



Figure 4 Flow chart and T-s diagram of the Stirling engine under finite 
heat reservoirs, (a) Simplified flow chart of the Stirling cycle under finite heat 
reservoirs; (b) T-s diagram of the Stirling cycle under finite heat reservoirs. 


are dumped into the environment at temperature T 0 . 

For such an engine, the heat absorbed from the hot 
stream and that released to the cold stream by the working 
fluid during a cycle are 

Qh = Cn^H-in -Tu-oJh = -r,)Ti, (26) 

Ql = C l (T l . oM -T l _Jt 2 = £ L C L _ mm (T 2 (27) 

where 7]and the 7) _ olli are the outlet temperatures of the 
hot and cold streams, respectively, and Cn-min and C L _ m j n are 
C H -n,in = min {C H , C g }, C L . min = min {C L , C g }. (28) 

where C H , C L and C g are the heat capacity flow rates of the 
hot stream, the cold stream and the working fluid through 
the heat exchangers, respectively. As the working fluid ex¬ 
periences the isothermal processes in the exchangers, C g can 
be considered as infinity. Therefore, Cn-min = Ch and CL-min = 
Cl- The effectivenesses of the exchangers on the hot and 
cold sides, e H and s L , can be calculated by 

^ H =l-exp[-(CA) H /C H . min ], (29) 

G. = 1 - exp [- (CA) L /C Lmin ], (30) 

where (t/A) H and (UA) t are the thermal conductance of the 
heat exchangers on hot and cold sides, respectively. 

Combining eqs. (5), (6), (26) and (27), the times of pro¬ 
cess a-b and process c-d can be calculated, 

1 ^.^(^-7;) 


iiRT 2 In A + »c v (1 - n m X?i ~ T 2 ) 
£ L C L-mJ T 2- T L-J 


The times of regeneration processes a-d and b-c, r 3 and 14 , 
can also be obtained by eqs. ( 11 ) and ( 12 ), and then the cy¬ 
cle time can be obtained by eq. (13). 

We use a similar isochoric temperature ratio exponent r 
with which the working temperatures 7) and 77 and the inlet 
temperatures Tu. m and satisfy 

^i=(7l in /7H.J- (33) 


According to the analyses above, for arbitrary r, all the 
parameters of the system can be obtained based on eqs. (5), 
( 6 ), and (26)-(33) with different 7). Then, the output power 
of the cycle is 

W =(Q h -Q L )/T = [nRlnA(T I -T 2 )\/t. (34) 

The entransy and entropy analyses are carried out below. 
Dividing the EL of each cycle by the cycle time, the EL rate 
can be obtained as 


G. oss =^c H cr„ 2 . m -T 0 2 )+^C l ( 7^ -T 0 2 ) 

4Q (r H _ out - T 0 )r, + C H (7^ - r 0 )(r - Tj) 

+ C L (T L . out -T 0 )T 2 +C L (T L _ m -T 0 )(t-t 2 )]—, (35) 

T 
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where the EL induced by dumping the used streams into the 
environment is considered [54], The first two terms on the 
right hand side are the decrements of the entransy flow rate 
of the hot and cold streams from the inlet to the environ¬ 
ment, and the third term is the increment of entransy flow 
rate of the environment. In one cycle time interval, r, the 
heat transfer between the hot stream and the working fluid 
only takes zi, the heat transfer between the cold stream and 
the working fluid only takes r 2 , and the streams dump to the 
environment directly in the rest time. Eq. (35) can be sim¬ 
plified into 

G loss = \c H (Tl in -r 0 2 )+ic L (r L 2 . in - T 0 2 ) 

-C H (T H _ m -T 0 )T 0 -C L (T L _ in -T Q )T 0 
[C H ( 7 H-i, - Th-ou,)T + C L (T L _ m - r L _ out )r 2 ] T 0 

T 

= 7 C H ( 7 H-in -T 0 f + \c L (T L . in -T 0 ) 2 +WT 0 . ,M ' ! 

It can be seen that the EL rate would increase monotonously 
with increasing output power when the heat capacity flow 
rates, the stream inlet temperatures and the ambient temper¬ 
ature are prescribed. 

Dividing the total EG of each cycle by the cycle time, the 
EG rate is 

S = C H ln—+ Cj In—2— 

gen n rji L rji 

*H-in *L-in 

+ [C H ( r H-o* - T oX+C H (T H _ m -T^z-z,) 

+ C L (^L-ou. - T 0 )G + C L ( r L-i„ - T 0 )(? - T 2 )]/(^0 ) > (37) 


where EG induced by dumping the used streams into the 
environment is considered [36]. The first two terms on the 
right hand side are the entropy variation rates of the hot and 
cold streams from the inlet to the environment, respectively, 
and the third term is the entropy variation rate of the envi¬ 
ronment. Eq. (37) can be simplified into 

^ g en = C n In + C]_ In z~~ 

-*H-in *L-in 

, C H (T H _ m -T 0 ) + C L (T L _ m -T 0 ) 

To 

, O^H-out -^H-in)G + Q.( 7 l-ou, ~ T L-Jh 
*To 

= C H ln^-+C L ln^- 
T\i-in T L -in 

,C H (T HJm -T 0 ) + C L (T LAn -T 0 ) W 

+- j - y (38) 

■'o 0 

It can be seen that EG also decreases monotonously with 
increasing output power when the heat capacity flow rates, 
the inlet temperatures of the streams and the ambient tem¬ 


perature are prescribed. 

Similarly, dividing the total ED of a cycle by the cycle 
time, the ED rate is 

c dis = \c a (Tl in -r 0 2 )+lc L (r L 2 in -r-) 

-[(QT, -Q 2 TJ/z] 

- [C H (7h. ou , - T 0 )z , + C H (r H . in - T 0 )(z - z,) 

+ Q (T L . 0Ut - T 0 )z 2 + C L (T LAn - T 0 )(z - z 2 )] ^, (39) 

r 

where ED induced by dumping the used streams into the 
environment is also considered. Eq. (39) can also be simpli¬ 
fied into 

G dis = \c^(Ty -7; 2 ) + Ic L (7’ L 2 in -T 0 2 ) 
-C H (T H _ m -T 0 )T 0 -C L (T LAn -T Q )T 0 
-[(QA-m/z] 

+ {K(T H _ in - T H out )z t +C L (T L m -T l _ 0 Jz 2 ]T q }/z 
= \c H (T H _ m -T 0 y+±C L (T L _ m -T 0 ) 2 
-[(Q 1 T 1 -Q 2 T 2 )/z] + WT q . (40) 

There is no specific relationship between the variations 
of the ED rate and the output power even when all the inlet 
parameters of the streams are prescribed. 

Numerical examples are given in the following paragraph 
to illustrate the applicability of EL, EG and ED to the opti¬ 
mization of the system with the objective of the maximum 
power output. 

The case in which 7) is optimized under the prescribed r 
is discussed first. Assume that r = 0.6, r H .i n = 800 K, 7j _ in = 
290 K, T 0 = 300 K, C H = 40 W K _1 , C L = 100 W K _1 , (CM) H 
= 50 W K -1 , (UA) L = 50 W K _1 , rj m = 0.8, c v = 20 J (mol K) _1 , 
R = 8.314 J (mol K) _1 , X = 3, n = 0.5 mol and Af, = M 2 = 
8x10 4 K s - 1 . According to the above parameters and equa¬ 
tions, the output power, the EL rate, the EG rate and the ED 
rate with different T) are calculated as shown in Figure 5. It 
can be seen that the output power and the EL rate both in¬ 
crease first and then decrease with increasing 7), while the 
EG rate and the ED rate both decrease first and then in¬ 
crease. When 7) = 660 K, the output power and the EL rate 
reach their maximum values, and the EG rate reaches its 
minimum value. However, the ED rate reaches its minimum 
value when 7) = 679 K. Thus, the maximum EL rate and the 
minimum EG rate both correspond to the maximum output 
power for the present case, while the extremum ED rate 
does not. Therefore, the maximum EL rate and the mini¬ 
mum EG rate both can be used in the optimization of this 
case, which is unanimous with the analyses above. 

For arbitrary r, the optimal working temperature 7] that 
leads to the maximum output power, the maximum EL rate 
and the minimum EG rate can be found. Therefore, a case 
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Figure 5 Variations of the output power, the EL rate, the EG rate and the 
ED rate with T\. 

with the variation of r is discussed next. For each r, the 
maximum output power, the maximum EL rate and the 
minimum EG rate at the optimal T) are calculated and 
shown in Figure 6. The maximum output power and the 
maximum EL rate both increase first and then decrease with 
increasing r, while the minimum EG rate decreases first and 
then increases. When r = 0.45, the output power and the EL 
rate reach their maximum values, and the EG rate reaches 
its minimum value. Therefore, the maximum EL rate and 
the minimum EG rate both can be used for the optimization 
of r to obtain the maximum output power. 

Moreover, the EGN of this problem can be calculated by 
eq. (23), where C mm = minjCn, C g , Cl}. According to eq. 
(24), the MEGN is 

= S gm T 0 /[C H (T H _ m —T 0 ) + C L (r L . in - r 0 )]. (41) 

Eqs. (23) and (41) reveal that the EGN and the MEGN both 
have the same variation tendencies with the EG rate when 
the inlet parameters (temperature, heat capacity flow) of the 
streams and the ambient temperature are prescribed, which 
demonstrates that both the EGN and the MEGN can be used 
for the power output optimization under the conditions 
above. 



Figure 6 Variations of the maximum output power, the EL rate and the 
minimum EG rate with r. 
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Another case in which the heat capacity flow rate of the 
hot stream Ch is not prescribed is considered. In this case, 
the total heat flow rate into the system is not constant. The 
working temperature is fixed at I) = 656 K and other pa¬ 
rameters are the same as the above numerical example. The 
output power, the EL rate, the EG rate and the ED rate are 
calculated and shown in Figure 7. All the parameters in 
Figure 7 increase monotonously with increasing Ch, which 
means that larger EL rate and larger ED rate both corre¬ 
spond to larger output power, while smaller EG rate does 
not. Therefore, the maximum EL rate and the extremum ED 
rate are suitable for the optimization of this case, but the 
minimum EG rate is not. Furthermore, the variations of 
EGN and MEGN with Ch are shown in Figure 8. It can be 
seen that EGN and MEGN both decrease first and then in¬ 
crease with increasing Ch- As the output power increases 
monotonously with increasing Ch, neither of EGN and 
MEGN fits the output power optimization when Ch is not 
prescribed. 

Actually, EG is a physical quantity that measures the ir¬ 
reversibility of processes. EG can be affected by two im¬ 
portant factors in heat-work conversion. One is the temper¬ 
ature difference of the heat transfer, and the other is the total 



Figure 7 Variations of the output power, the EL rate, the EG rate and the 
ED rate with Ch. 



Figure 8 Variations of the output power, the EGN and the MEGN with 

Ch- 
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heat transferred. When the inlet parameters of the hot 
streams are prescribed, the total heat that can be used for 
doing work is prescribed. Under this condition, smaller EG 
rate means less loss of the ability of doing work, which 
leads to larger output power. When the inlet parameters are 
not prescribed, the total heat that can be used for doing 
work is not prescribed. The output power of the system is 
not only related to the temperature difference of the heat 
transfer, but also to the total heat input. With the increase in 
the total heat input, the output power increases, and the EG 
rate increases too. Consequently, the minimum EG rate is 
not a suitable criterion for the output power optimization in 
such a condition. Overall, EG does not describe the output 
power of the system directly, but indirectly by the loss of 
the ability of doing work. Just due to this indirectness, limi¬ 
tations are encountered when EG is used for optimizing the 
thermal processes. 

5 Conclusions 

The entransy theory is combined with the finite time ther¬ 
modynamic theory to analyze the Stirling engine cycles in 
this paper. The applicability of the entransy and entropy 
theories to the optimizations of the cycle with the objective 
of the maximum output power is investigated. 

(1) For the cycle working under the infinite heat reser¬ 
voirs whose temperatures are prescribed, the maximum EL 
rate can be used for optimizing the power output. The min¬ 
imum EG rate is not related to the maximum power output; 
instead, the maximum EG rate corresponds to the maximum 
output power. The extremum ED rate has no direct relation 
to the power output and is not appropriate for the optimiza¬ 
tion. EGN cannot be defined, and MEGN also has no direct 
relation to the output power because it is a constant when r 
is prescribed. 

(2) For the cycle working under the finite heat reservoirs 
provided by the streams, the maximum EL rate, the mini¬ 
mum EG rate, the minimum EGN and the minimum MEGN 
all can be used for the output power optimization when the 
stream inlet temperatures and the heat capacity flow rates 
are prescribed except the extremum ED rate. When the heat 
capacity flow of the hot stream is not prescribed, the mini¬ 
mum EG rate does not corresponds to the maximum output 
power, but the maximum EL still does. 

(3) The concept of entransy loss considers the effects of 
the heat transfer processes and the work doing processes on 
the variation of the entransy, i.e., ED and work entransy. 
For the cases discussed in this paper, the increase in heat 
transfer rate leads to the increase in output power and the 
increase in the ED rate, and finally results in the increase in 
the EL rate. Similarly, the increase of the output work leads 
to the increase in work entransy flow rate, and then the EL 
rate also increases. Thus, a better consistency is discovered 
when the concept of EL is used for the output power opti¬ 


mization of the Stirling cycle. For the concept of EG, larger 
output power corresponds to the smaller EG rate when the 
total heat input from the hot stream is prescribed. When the 
total heat input from the hot stream increases under the pre¬ 
scribed temperature, the output power increases and the EG 
rate increases, too. The relationship between the variations 
of the output power and the EG rate is affected by the pre¬ 
conditions, which makes the applicability of EG, EGN and 
MEGN conditional. The concept of ED is not always ap¬ 
propriate for the heat-work conversion optimization because 
it does not consider the effect of the work doing processes 
on the variation of the entransy. 
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